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It is well known that a random multiplicative process with 
weak additive noise generates a power-law probability distri- 
bution. It has recently been recognized that this process ex- 
hibits another type of power law: the moment of the stochas- 
tic variable scales as a function of the additive noise strength. 
We clarify the mechanism for this power-law behavior of mo- 
ments by treating a simple Langevin-type model both approx- 
imately and exactly, and argue this mechanism is universal. 
We also discuss the relevance of our findings to noisy on-off 
intermittency and to singular spatio-temporal chaos recently 
observed in systems of non-locally coupled elements. 
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I. INTRODUCTION 

Power laws are observed in a wide variety of natural 
phenomena and mathematical models. Some examples 
are the critical behavior near second order phase tran- 
sitions, Kolmogorov's law of fully developed turbulence, 
size distribution of avalanches in models of self-organized 
criticality, Gutenberg-Richter's law of earthquakes, dis- 
tribution of price fluctuations in economic activities, and 
Zipf's law in linguistics. Clarifying the mechanisms for 
the emergence of these power laws has long been a sub- 
ject of many challenges. 

The random multiplicative process (RMP) is a well- 
known mechanism leading to power-law behavior. It is a 
stochastic process where the stochastic variable is driven 
by a multiplicative noise. It has been extensively used as 
a model for a variety of systems such as on-off intermit- 
tency lasers [Q, economic activity H||, variation 
of biological populations in fluctuating environment JTo[ , 
and passive scalar field advected by fluid . 

In real systems, the stochastic variable may often be 
driven not only by the multiplicative noise, but also by 
some weak additive noise. This weak additive noise be- 
comes important when the amplitude of the stochastic 
variable x takes small values, and introduces an effec- 
tive lower bound of x. Actually, this lower bound may 
be crucial, because it guarantees the existence of a sta- 
tionary probability distribution function (PDF). Further- 
more, the PDF here has a power-law form over a wide 
range of a; 0|g§§ll]|| . 

For example, Venkataramani et al. [|5| introduced a 
Langevin equation with multiplicative and additive noise 
terms as a model for noisy on-off intermittency. They 
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obtained a stationary PDF with a power-law tail. The 
same form of Langevin equation was treated by Takayasu 
et al. ||] as a model for economic activity, and they also 
showed that the PDF obeys a power law. A similar model 
was introduced by Levy et al. H: it describes a discrete 
stochastic process driven by a multiplicative noise. They 
introduced a lower bound to the stochastic variable ex- 
plicitly, and showed again that the PDF obeys a power 
law. Venkataramani et al. and Takayasu et al. treated 
the additive noise explicitly, while the lower bound intro- 
duced by Levy et al. plays a role similar to the additive 
noise. In this paper, we are concerned with this type of 
stochastic processes. 

Recently, another type of asymptotic power law was 
found in the above type of stochastic processes. In previ- 
ous papers 0,[D|], we introduced a stochastic process in 
order to explain the power law displayed by the spatial 
correlation function C(r), i.e., C(r) ~ Co—Cir a for small 
enough r, observed in the spatio-temporal chaotic regime 
of systems with non-locally coupled elements. Our ex- 
planation was based on a RMP with weak additive noise 
such as described above. Note, however, that the power- 
law correlation here is not a direct result of the power-law 
tail of the PDF itself, but it is a result of the asymptotic 
power law of moments < x q > of the stochastic variable 
x as a function of the strength s of the additive noise, 
i.e., < x q >~ Go + G\s H ( q \ This gives another mech- 
anism leading to power-law behavior in such stochastic 
processes. 

The goal of this paper is to clarify this mechanism for 
the emergence of the power law of moments with respect 
to the strength of the additive noise. We achieve this 
by using a simple Langevin-type model, and argue that 
the mechanism proposed is a universal one in generating 
various power laws. 

The outline of this paper is as follows: In Sec. II we 
introduce the model to be studied, and display its typical 
behavior. In Sec. Ill we treat the model approximately 
in order to outline the mechanism for the emergence of 
the power law of moments, and then exactly in Sec. IV. 
In Sec. V we discuss the robustness of the power law 
with regard to boundary conditions and nature of the 
noise. We also show some results obtained by numerical 
calculations with colored noises. In Sec. VI we discuss 
an application of our theory to noisy on-off intermittency. 
As an example, we show a result obtained by a numerical 
calculation of coupled chaotic elements. Furthermore, 
we discuss the relation of the power law of moments to 
the power-law spatial correlations observed in systems of 
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non-locally coupled chaotic elements. We summarize our 
results in Sec. VII. 



II. ANALYTICAL MODEL 

A. Langevin equation 

As a model for a RMP with weak additive noise, we 
employ a Langevin equation 

*& = \(t)x{t)+ V (t), (1) 

where x(t) is a stochastic variable, X(t) a multiplicative 
noise, and rj(t) an additive noise. We assume both X(t) 
and rj(t) to be Gaussian-white, and their average and 
variance to be given by 

< \(t) >= A , < [A(t) - A ] [A(f ) - A ] >= 2D x 8{t - t'), 
< r){t) >= 0, < T){t) rj(t') >= 2D v 5(t - t'). 

(2) 

We further assume D v <C D\, namely, the additive noise 
is sufficiently weaker than the multiplicative noise. 

This simple Langevin equation ([!]) has been widely 
used in many studies of various systems PM^O] . The 
physical meaning of x(t), X(t), and n(t) may be different 
depending on the specific system under consideration. 
For example, in the case of lasers, x(t) gives the number 
of photons, X(t) its fluctuating amplification rate, and 
•q(t) the noise due to random spontaneous emissions of 
atoms. When we are working with noisy on-off inter- 
mittency, x(t) gives the measure of a distance from the 
invariant manifold, X(t) the instantaneous vertical Lya- 
punov exponent, and rj(t) the noise due to a parameter 
mismatch or some other cause. In the context of eco- 
nomic models, x(t) represents the wealth, X(t) the rate 
of change of the wealth, and ij(t) some external noise of 
various sources. 



B. Boundary conditions 



(b) Upper bound To be realistic, when x(t) takes too 
large values, it should be saturated by some effect such as 
nonlinearity. We simply introduce this effect as boundary 
conditions, specifically reflective walls at x = ±1. 

With these upper and lower bounds provided by an 
additive noise and reflective walls, the Langevin equa- 
tion (|l|) admits a statistically stationary state. In Fig. [I], 
we show a typical time evolution of x(t) governed by the 
Langevin equation (|l|) for slightly negative Ao, where the 
reflective walls are placed at x = ±1. In spite of the 
negative average expansion rate Ao, x(t) does not simply 
decay but exhibits intermittent bursts. The generation 
of bursts may be interpreted as follows. Due to the weak 
additive noise, x(t) may generally have small but finite 
values. If positive X(t) happens to persist over some pe- 
riod, x(t) will be amplified exponentially and attain large 
values, which is nothing but bursts. Of course, x(t) may 
eventually decay to the noise level because of negative Ao . 
The chance of bursts will increase with the additive noise 
strength. Why this leads to a power-law dependence of 
moments of x(t) on the additive noise strength can be 
understood from the argument below. 

The intermittency described above has the same sta- 
tistical nature as the noisy on-off intermittency. Actu- 
ally, one may consider the noisy on-off intermittency as 
a stochastic process of the type described above. 
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In order to obtain a statistically stationary state from 
the Langevin equation (|l|), we generally need upper and 
lower bounds of x. In our model, the weak additive noise 
may act as an effective lower bound. 

(a) Lower bound Without the additive noise, x(t) 
tends to when the average expansion rate Ao is neg- 
ative. The additive noise introduces an effective lower 
bound of x(t), which keeps x(t) away from the zero value 
even if Ao < 0. In this paper, we treat the additive noise 
explicitly, while in some other studies it is replaced by 
a reflective wall (infinitely high barrier) placed at some 
small x. 



FIG. 1. Typical time evolution of the amplitude x(t) of 
the Langevin equation (1) for Ao = —0.5, D\ = 0.5, 
D v = 0.00005. 
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III. APPROXIMATE TREATMENT 



(9) 



In order to give an outline of the mechanism underlying 
the emergence of the power law of moments, we first treat 
the Langevin equation ([!]) approximately. 



A. Fokker-Planck equation 



We introduce a characteristic amplitude 



where the flux j2{x, t) is now defined by 
j 2 (x,t) = -D- n ^-P 2 {x,t). 



(10) 



The boundary conditions to be imposed here are the 
continuity of Pi and P 2 , and also j\ and j 2 , each at 
M =s (y = logs). 




(0 < s < 1) 



(3) 



B. Stationary PDF with a power-law tail 



which results from the balance between the fluctuation 
due to the multiplicative noise < (Aa;) 2 >~ D\x 2 and the 
fluctuation due to the additive noise < rj 2 >~ D v . We 
divide the range of x into two parts, < \x\ < s and 
s < |ar| < 1, and ignore in each region one of the noise 
sources which is less dominant there. Since the system 
is statistically symmetric with respect to the transforma- 
tion x — » —x, we consider only the absolute value |x| 
hereafter. 

(a) s < \x\ < 1 In this region, we ignore the effect of 
the additive noise and consider a Langevin equation 



dx(t) 
dt 



\{t)x(t). 



(4) 



By introducing a new variable y(t) = log|x(i)|, Eq. 
is rewritten as 



dy(t) 

dt 



A(i), 



(5) 



which gives a diffusion process with mean drift. 

Let Pi(x,t) denote the PDF of x{t) and P x {y,t) the 
corresponding PDF of y(t). The Fokker-Planck equation 
corresponding to Eq. (0) takes the form 



where the flux j\{y, t) is given by 



(6) 



j 1 {y,t) = \ Q Pi{y,t)-D x —Pi{y,t), (7) 

dy 

Setting reflective walls at x = ±1 is equivalent to assum- 
ing a no-flux boundary condition for j\ at y = 0, i.e., 

ii(» = o,t) = o. 

(b) < \x\ < s In this region, we ignore the effect 
of the multiplicative noise, and this leads to a Langevin 
equation 



dx(t) 
dt 



= ri(t). 



(8) 



Let P 2 (x,t) be the PDF of x(t). P 2 {x,t) obeys the 
Fokker-Planck equation 



We calculate here the stationary solution of the Fokker- 
Planck equation. 

(a) s < \x\ < 1 Stationarity condition dP\/dt = 
gives dji/dy = 0, i.e., ji(y) = const., and the no-flux 
boundary condition ji(y — 0) = gives ji(y) = 0. 
Therefore, the stationary solution P\ (y) satisfies 



= XoPi(y) - D x -^Pi(y). 



This can be solved as 



/» (//) = T'rxp ( ^-y 



(11) 



(12) 



where C is a normalization constant. In terms of the 
original variable x, we obtain 



P 1 (x)=P 1 (y) 



dy 
dx 



CjH-exp ( ^-l g(|x|) 
\x\ \Dx 

C\x\~ \ 



(13) 



Thus, the PDF obeys a power law in this region. The 
exponent of this power law is determined by the ratio of 
Ao to D\, i.e., by the basic statistical characteristics of 
the multiplier A(i), and does not depend on the nature of 
the additive noise. We denote this ratio as (3 hereafter: 



(14) 



(b) < \x\ < s A general form of the stationary 
solution is given by P 2 (x) — Ax + B, where A and B are 
constants. Continuity of the flux and the PDF at |x| = s, 
i.e., j 2 {s) — ji(s) = and P 2 (s) = Pi(s) gives A = and 
B = Pi(s). Therefore, ^2(2;) takes a constant value: 



P 2 {x) = Pi(s) = Cs^ 



(15) 



Finally, the approximate stationary PDF is obtained 



as 
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CV 3 - 1 (0 < \x\ < s) 
P(x) = { Cx?- 1 (s < |x| < 1) • 
(|a:| >1) 
The normalization constant G is determined from 



P(x)dx = 2 P(x)dx = 1, 
Jo 



(16) 



(17) 



and calculated as 



C 



x - x dx 



2 s? + 







(18) 



Thus, the PDF consists of three parts, i.e., a constant 
part near the origin where x(t) is dominated by a normal 
diffusion process, a power-law tail where x(t) is domi- 
nated by a RMP, and a vanishing part. The boundary 
between the constant part and the power-law tail is lo- 
cated at | a; | = s, which is proportional to the additive 
noise strength JD^. We show this approximate PDF 
@ in Fig. | 
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FIG. 2. PDFs P(x) vs. x, where the approximate and 
exact curves are compared. The parameters are the same as 
in Fig. 1. 



C. Moments 

The q-th moment < x q > of |arj in this stationary state 
(O) is calculated as 



<x q > = 



J \x\ q P{x)dx = 2 J x q P{x)da 



= 2C(J x q s^ 1 dx + J x q x^~ x dx 

p ;+(f±g-i)A 

+ q l + (/3-l)s 
We can write the above in the form 

1 + ClqsP+l 



< X q >- 



1 + olqsI 3 



where c q and a q are given by 



(3+1 
l + o 



1. 



(19) 



(20) 



(21) 



Note that the exponent (3 of the PDF now appears as the 
power of s. As we explain below, the form of Eq. ( |20| ) 
is all we need for the emergence of the asymptotic power 
law of moments. 



D. Asymptotic forms of the moments 

We investigate the asymptotic forms of the moment 
< x q > in the limit of small additive noise s. We consider 
only the practically interesting case of positive q. 

(a) P > By expanding the denominator of < x q > 
and taking the lowest order in s, we obtain 



< x- 1 >: 



(l - a 



(22) 



(b) (3 < Ignoring 1 in the denominator of < x q >, 
we obtain 



< x H >: 



2l( s W +aqS q ) 



(23) 



Which of \(3 1 and q is smaller determines which of the two 
terms on the right-hand side of Eq. (|2^) dominates. We 
thus obtain 



< x* >: 



tl s \P\ 
a Q 



Cqa,q gg 
a 



><?) 



(24) 



These results show that the moment < x q > ap- 
proaches a simple power-law form as a function of the 
position s of the boundary, or the strength yj D v of the 
additive noise: 



G + G lS H ^\ 



(25) 



where Go and G\ are constants. Go vanishes when j3 < 
(Ao < 0), while it takes a finite value when f3 > (Aq > 
0). Thus, we have obtained the power law of moments. 
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E. Exponents 



IV. EXACT TREATMENT 



The exponent H(q) of the moment < x q > is deter- 
mined by /3, namely, the ratio of Ao to D x . From Eqs. 
( p2| ) and (|24|), H{q) varies with g as follows: 

(a) (3 < 



H(q) = 



q (0< g < 



(g) = /3. 



(26) 



(27) 



We notice that H(q) = q when < q < |/3|, but i?(g) = 
|/3| when |/3| < g or j3 > without dependence on q. 



F. Other asymptotic regimes 

When j3 ~ or |/3| ~ g, there exist other asymptotic 
regimes where the asymptotic form of the moment in the 
s — > limit is not a power law. 

(a) (3 ~ Consider a parameter region where /3 ~ 
and |/3 logs | <C 1. The denominator of < x q > can be 
expanded as 

1 + a s /3 = 1 + a exp(/3 log s) 

= l + a + a /31ogs + C(|/3 log s| 2 ) . (28) 

Using ao = f3 — 1 and | logs| 3> 1, the ao/31ogs term is 
found to be dominant and we obtain 



\a /3 log s\ 



(29) 



Thus < x q > diverges logarithmically as 1/| logs|. 

(b) |/3 1 ~ q Consider a parameter region where /3 < 
and |/3| ~ g. We further assume |(g — |/3|)logs| <C 1. 
< x q > is then given by 



< x* >c 



2i( s \p\ + a s i) = 2i s m (i + a qS q 

ao V / ao V 



-|/3| 



(30) 



Expanding the right-hand side, we obtain 



I +a q s q 1/51 = l + a 9 cxp((g- |/3|)logs) 
= l + a g + a,(g- |/3|)logs 

+ 0(|(g-|/3|)log S | 2 ) 



By using a q — f^ — 1 and | logs| » 1, the a 9 (g- 
term is found to be dominant, and we obtain 

^ n da 



< x q >c 



ao 



|(g-|/3|) S l^log 



s . 



(31) 

! logs 

(32) 



Thus < x q > diverges as js'' 3 ' logs| 



Next, we treat the effect of the additive noise without 
approximation. The argument below is in parallel with 
the previous one, and their results agree qualitatively, 
giving the same values of the exponents. How to calcu- 
late the PDF follows the argument by Venkataramani et 

al. a. 



A. Fokker-Planck equation 

The Fokker-Planck equation corresponding to the 
Langevin equation ([!]) is given by 



^p(*,t) = -A i(a , t)> 



(33) 



where P(x, t) is the PDF of x(t), and the flux j(x, t) takes 
the form 

j(x,t) = (A + D x )xP(x,t) - -^[(D x x 2 + D n )P(x,t)]. 

(34) 

Reflective walls at x = ±1 are equivalent to imposing 
no-flux boundary conditions at x = ±1, i.e., 



j(x=±l,t)=0. 



(35) 



B. Stationary PDF with a power-law tail 

We calculate the stationary solution P(x) of the 
Fokker-Planck equation (|33|). Stationarity condition 
dP(x, t)/dt = gives dj(x, t)/dx — 0, i.e., j(x) = const., 
and no-flux boundary conditions j(x = ±1) = give 
j(x) = 0. Therefore, P{x) obeys 

(Ao + D x )xP{x) - -^[(D x x 2 + D n )P{x)] = 0. (36) 



Solving this, we obtain 



P(x) = C{D x x 2 + D rt 



(37) 



as the stationary PDF, where C is a normalization con- 
stant to be determined from 



P{x)dx = 2 / P(x)dx = 1. 



-1 

If we use the integral formula 

r 1 

„2\a 



1 3 



(l + cx 2 ) a dx = 2 F 1 (-a,-,-;-c), 
o l i 



(38) 



(39) 



where iF\{a, b, c; z) is the hypergeometric function, C can 
be expressed as 
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c = 



2d; 



2^1 (- 



Aq 
2D, 



1 1 3._-Da, 
2 ' 2 ' 2 ' D„ ' 



(40) 



We define /? as the ratio of the average expansion rate Ao 
to its fluctuation D\ as in the previous calculation, and 
a as the exponent of P(x), i.e., 



a 



(3-1 



Aq 
2D. 



(41) 



numerically by a direct simulation of the Langevin equa- 
tion (0). Figure || illustrates PDFs for different values 
of Ao with fixed D v , while those in Fig. [| are for dif- 
ferent values of D v with fixed Ao- Each PDF takes a 
constant value near the origin, whereas it obeys a power 
law otherwise. Their exponents vary with Ao, and the 
crossover position moves to the right with the increase of 
the strength 



'D v of the additive noise. 



Further, we define s as the ratio of the strength y/D v of 
the additive noise to the strength \JD\ of the multiplica- 
tive noise: 



(42) 




Finally, the stationary PDF is expressed as 
(l + f^) Q 



P(x) 



22Fi(—a, j, §; 







(N<i) 



X > 1 



(43) 



This stationary PDF approaches a constant as x — > 0, 
and when s 2 <C 1, namely, when the additive noise is 
sufficiently weaker than the multiplicative noise, the PDF 
approaches a power law as x — ► ±1: 



P(x) 



const, (x — > 0) 
\x\ 2a (x^±l) 



(44) 
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FIG. 3. PDFs P(x) vs. x obtained numerically(o) and 
theoretically ( — ) for different values of Ao- D\ — 0.5 and 
D v = 0.00005 are fixed. 



The exponent of this power law is given by 



Ao 



2a = 0-1= 



(45) 



Thus, we obtain qualitatively the same PDF as the pre- 
vious approximate result ([lq), in particular, a power law 
with the same exponent. The crossover point between the 
constant region and the power-law region of the PDF is 
found from the balance of the two terms in the numerator 
of P(x): 



s 2 ' 



The crossover thus occurs near 



(46) 



100 r 



0.01 



0.001 




(47) 



which is exactly the point at which we divided the domain 
of x in the previous approximate treatment. 

We show the exact PDF @ in Fig. |. The approxi- 
mate PDF ([l6]) reproduces the main features of the exact 
one well. In Fig. | and^, we display two graphs of PDFs, 
one obtained theoretically in (E3|) and the other obtained 



FIG. 4. PDFs P(x) vs. x obtained numerically (o) and 
theoretically ( — ) for different values of D v . Ao = —0.5 and 
D\ = 0.5 are fixed. 
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C. Moments 

The q-th moment < x q > with respect to the station- 
ary PDF (43) is calculated as 



< x q > 



\x\ q P{x)dx = 2 [ x q P(x)dx 
1 Jo 

1 TP ( ™ 1 +1 3+q. 1 \ 



1 + q (-a, 
where we used the integral formula 



(48) 



z 9 (l + ca; 2 ) a dx = 
o 1 + 5 



1 1+g 3+g 

2-Fi(-a, — ^— , — ^— ; -c). 



(49) 

In Fig. ||, we show the moments < x q > obtained theo- 
retically in ( |48"| ) and compare them with those obtained 
numerically. 



0.1 : 



0.01 : 




< x q >: 



1 ri (-a,i±2,5±2) l + r 3 (-a, , ■ 



1±3 3 +9 ) s /3+g 



1 + q IM-ai, 



1 3> 

2 ' 2 - 1 



i + r 3 (-a ^ 



2' 2' 

(51) 

where T\, and r 3 are defined in terms of the gamma 
function T(a) as 

r(c)r(6- a ) , r(c)r( -h) 

ri(a,6,c) = r, r 2 (a,&,c; - 



and 



r(&)r(e-a) 



r 3 (a,6,c) 



r 2 (a,b, c) 
ri(a,6, c)' 



r(o)r(c-6)' 

(52) 
(53) 



Notice that here again we obtain the form already ob- 
tained in the previous approximate calculation: 



< x H >= c, 



1 + a qS P +q 



1 + a s 







(54) 



However, the expressions for c q and a q are different, and 
now given by 



1 rj(-a,i±g,2± g) 
1 + Q T^-a,!,!) 



. ■ j ; r / 1 + g 3 + g 

= r 3 (-a, — r— , — — ). 



(55) 



Using this form, and from exactly the same reasoning as 
before, we can show that < x q > asymptotically obeys a 
power law as s — ► 0: 



<x q >~ G + G lS H{q) . 



(56) 



In Fig. [3], we show the moments < x q > for small s ob- 
tained both theoretically and numerically. Each moment 
shows power-law dependence on the strength s of the ad- 
ditive noise. 



FIG. 5. Moments < x q > vs. the strength s of the additive 
noise obtained numerically (symbols) and theoretically(lines) 
for several values of q. The parameters are Ao = —0.5 and 
D x = 0.5. 



D. Asymptotic forms of the moments 

We investigate the asymptotic forms of the moment 
< x q > in the limit of small additive noise, s — > 0. As 
before, we consider only the case q > 0. 

Using the asymptotic form of the hypergeometric func- 
tion, i.e., 

2 F\(a, 6, c; z) ~ Ti(a, b. c)(-z)~ a + T 2 (a, 6, c)(-z)~ 6 

(z -> 00), (50) 
we can write the asymptotic form of < a; 9 > as 



E. Exponents and other asymptotic regimes 

The only difference between the above exact result and 
the previous approximate result is in some coefficients in- 
volved. Since the exponents are unchanged, the behavior 
of H(q) is exactly the same as the previous result. There 
is also no difference in that there exist other asymptotic 
regimes near /3 = or \j3\ = q, which is clear if we notice 
a q = -1+ {/3) and | log s| > 1. 

In Fig. [7], the exponent H(q) versus Ao obtained the- 
oretically in (|26| ) and (|27]) are given in comparison with 
those obtained numerically. Each H(q) — Xq curve is 
composed of two parts, i.e., a part where H(q) varies 
in proportion to | Aq | and that where H(q) saturates to 
a constant. We estimated the exponents numerically by 
assuming a power law even in the above-mentioned non- 
power-law asymptotic regimes. Therefore, the estimated 
values there are naturally different from those expected 
theoretically. 
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V. ROBUSTNESS OF THE POWER LAW 
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FIG. 6. Moments < x q > vs. the strength s of the additive 
noise obtained numerically (symbols) and theoretically(lines) ; 
a blowup of Fig. 5 for a small s region. 



Although we have treated only the Langevin equation 
(|l|) up to this point, the type of power law discussed 
above appears in many other models. It is insensitive to 
the details of the model such as the boundary conditions 
imposed, discreteness or continuity of time, and the na- 
ture of noise terms. We thus discuss the robustness of 
the power law here. 

A. Boundary conditions 

We treated the effect of the additive noise explicitly 
both in the approximate and exact calculations. The 
crucial role of the additive noise in generating a power- 
law PDF is to save the stochastic variable from decaying 
completely by generating small fluctuations around the 
zero value where a normal diffusion process dominates. 
This is the reason why the usual approximation of re- 
placing the additive noise with an explicit lower bound 
of the variable works well. 

Although we assumed that the upper bound of the 
stochastic variable is simply given by the reflective walls, 
the result would not change essentially if we replace it 
with some nonlinearity as given by a — x 3 term, at least 
for not too large q. This is because the dominant con- 
tribution to the s-dependence of < x q > comes from the 
region of large P(x), i.e., that of small x and not from 
the large x region near the upper bound. 



B. Discrete models 




FIG. 7. Exponents H(q) vs. Ao obtained numeri- 
cally (symbols) and theoretically(lines) for q — 1,2,3. The 
value of D\ is fixed to 0.5. 



Power laws also appear in discrete-time models, and 
their origin is exactly the same as before. For example, 
in 1 13 1 we introduced a discrete time stochastic process 



i"n+i 



0{xl 



Vn 



(57) 



where n is the time step, A„ and rj n are noise. We ap- 
proximated the additive noise term and the nonlinear 
term as lower and upper reflective walls, and obtained a 
power-law PDF. Furthermore, we obtained a power-law 
dependence of the moments of x on the position of the 
lower bound, i.e., the strength of the additive noise. 



C. Nature of noise 

Of course, assuming a Gaussian-white noise is some- 
times inadequate for models of real systems. The power 
law of moments can also be seen in some models with 
colored noises. For example, in fl^| we introduced a 
stochastic process driven by a colored dichotomous noise 
with lower and upper reflective walls. We obtained the 
power law of moments with respect to the position of the 
lower wall as well. 
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D. Numerical examples 



In order to demonstrate the robustness of the power 
law with regard to the nature of noise terms and bound- 
ary conditions, we give a few numerical examples below. 

We numerically study a stochastic process 

^ = [\ + \(t)}x(t)+r 1 (t), (58) 

where X(t) is some colored noise, while r)(t) is a Gaussian- 
white noise as before. We use three different types 
of noise for A(t): (i) Gauss-Markov noise produced by 
an Ornstein-Uhlenbeck process, (ii) dichotomous noise 
which takes values 1 or —1 with equal transition prob- 
abilities, and (iii) chaotic noise produced by the Lorenz 
model. We normalize the average and variance of each 
noise to and 1 respectively, and use it for A(t). We ad- 
just Ao to a value where x(t) shows intermittency similar 
to Fig. § 

In Fig. [8[ we show second moments < x 2 > versus the 
strength s of the additive noise obtained for each type 
of noise. We can see a power-law dependence of the mo- 
ments < x 2 > on the strength s of the additive noise in 
the small s region. We also studied the case where the 
saturation of x is not due to the reflective walls but a 
nonlinearity —a; 3 , and found that a power law with the 
same exponent still holds. We also confirmed that the 
PDF has a power-law tail for each type of noise. 

Since the noise is colored, it would be difficult to pre- 
dict the values of the exponents of the moments from the 
previous theory. In order to achieve this, some sort of 
renormalization procedure, like the one done in || , must 
be invoked to give effective Ao and D\. It is beyond the 
scope of this paper. 

VI. SOME RELATED SYSTEMS 
A. Noisy on-off intermittency 

Since noisy on-off intermittency is a typical phe- 
nomenon with the mechanism of generating the power 
law of moments, we briefly discuss it here. 

On-off intermittency is observed where a chaotic at- 
tractor becomes marginally stable with respect to dis- 
turbances transversal to the invariant manifold in which 
the chaotic attractor is embedded. This type of insta- 
bility is called a blowout bifurcation. The system then 
alternates between two phases intermittently: a laminar 
phase where the system stays practically on the invariant 
manifold, and a burst phase where the deviation from the 
invariant manifold grows suddenly. The mechanism re- 
sponsible for the on-off intermittency is that the distance 
between the orbit and the invariant manifold is governed 
by a multiplicative process with a chaotically changing 
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of the additive noise calculated for dichotomous noise, 
Gauss-Markov noise, and Lorenz noise. The value of Ao is 
—0.5 for the dichotomous and Gauss-Markov noise, and —0.3 
for the Lorenz noise. Each line is shifted upwards or down- 
wards to avoid overlap. 

multiplier. Therefore, the corresponding suitable math- 
ematical model is a RMP, where the fast chaotic motion 
of the multiplier is considered as a multiplicative noise. 

Piatt et al. B investigated a situation where weak ad- 
ditive noise is also present in the system. They found that 
the intermittency, which was originally observed only in 
a narrow supercritical-side region of the blowout bifur- 
cation, can be observed in a wider region including the 
subcritical side of the blowout bifurcation. This is called 
noisy on-off intermittency. In order to explain this phe- 
nomenon, they proposed a RMP model which incorpo- 
rates the effect of the additive noise approximately by 
introducing a lower bound to the conventional model of 
on-off intermittency. Venkataramani et al. |5j and Cenys 
et al. H also analyzed similar models and explained some 
features of the noisy on-off intermittency. However, the 
power law of moments which we described above has not 
fully been investigated. Thus, we give an example below. 

B. Coupled chaotic elements 

Two identical chaotic elements coupled with each 
other Q is a typical system which shows on-off inter- 
mittency. A frequently used model is 

X x = F(X 1 )+k(X 2 -X 1 ) 

(59) 

X 2 = F(X 2 ) + k(X 1 - X 2 ) 

where X = F(X) gives the dynamics of the individual 
element, and fc the coupling strength. These two ele- 
ments synchronize with each other when k is larger than 
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a certain critical coupling strength k c . The difference 
x = X2 — X\ is driven multiplicatively by the chaotic 
motion of the elements as 



2kl 



Q(x 2 ), 



(60) 



where D means differentiation and I the unit matrix. 
Slightly below k c , x exhibits on-off intermittency. If wc 
further apply some weak additive noise, x exhibits in- 
termittent behavior above k c as well. If we vary the 
strength of the additive noise and measure the q-th mo- 
ment < x q > of a certain component of x, it is expected 
to behave like < x q >~ Go + G\s H ^ as far as s is suffi- 
ciently small. 

As an example, we numerically calculate a pair of 
Rossler oscillators coupled with each other, and with 
weak additive noise applied only to the first component 
of the state variables: 



±i(t) 



-yi + zi + k(x 2 - xi) + s£i(t) 



and 



yi{t) = x 1 +Q.iyi + k{y 2 -yi) 

z\{t) = 0.2 + x\Z\ — 5.7zi + k(z2 — Z\) 



(1 — 2), 



(61) 



(62) 



where £1,2 (i) are Gaussian-white noise of average and 
variance 1, and s controls their strength. We set the 
coupling strength k to the value slightly above k Cl where 
the system shows noisy on-off intermittency. Figure ^ 
shows second moments < x 2 > of the difference between 
the first components x = x 2 — x\ versus the noise strength 
s, obtained for some values of the coupling strength k. 
As expected, they show power-law dependence on s, and 
their exponents vary with k. 

We also observed the power law of moments in a sys- 
tem of coupled chaotic elements with a slight parameter 
mismatch M in place of a weak additive noise, where the 
parameter mismatch plays a role similar to the additive 
noise. The power law of moments is also observed in 
systems of coupled maps [H. 



C. Spatially distributed chaotic elements 

The power-law spatial correlation which we studied in 
the previous papers Jl2],[l3| turns out to be the same as 
the power law described above, with a proper physical in- 
terpretation of some variables involved. The systems we 
treated previously are the populations of spatially dis- 
tributed chaotic elements. The elements are driven by a 
field produced by non-local coupling, which is spatially 
long- waved and temporally irregular. 

If we consider the difference x(r) = X(r) — X(0) be- 
tween the amplitude X of two elements which are sepa- 
rated by a short distance r, its moments < x(r) q > are 
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FIG. 9. Second moments < x 2 > vs. the strength s of the 
additive noise of coupled Rossler oscillators obtained for sev- 
eral values of k. 



directly related to the spatial correlations of X. They 
are similar to the structure functions in the context of 
fluid turbulence or the height-height correlations in the 
context of fractal surface growth. The dynamics of x(r) 
is given by a RMP with a weak additive noise, where 
the multiplier is the local Lyapunov exponent fluctuat- 
ing randomly due to the chaotic motion of the elements, 
and the weak additive noise comes from the small differ- 
ence in the strength of the applied field between the two 
points under consideration. 

Since the strength s of the additive noise should be of 
the order of r due to the assumed smooth spatial vari- 
ation of the applied field, the power law of moments as 
a function of s is now interpreted as a power law of mo- 
ments of the amplitude difference function of the 
mutual distance r, i.e., 

< x(r) q >~ Go + G lS H ^ ~ Go + G 2 r m<i \ (63) 

where Go, Gi, and G 2 are constants. 

When the average Lyapunov exponent is negative, Go 
vanishes and the exponent H(q) is given by Eq. (|2q) . 
It shows a "bifractality" similar to that known for Burg- 
ers equation, which implies an underlying intermittent 
structure. 



VII. CONCLUSION 

As in the models for noisy on-off intermittency and 
economic activity, a stochastic process driven by mul- 
tiplicative and weak additive noise shows a power-law 
PDF. The PDF consists of a constant part and a power- 
law part, and their boundary moves with the strength s 
of the additive noise in such a way that its distance from 
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the origin is proportional to s. This systematic depen- 
dence on s causes the power law of the moments with 
respect to the strength of the additive noise. 

In order to study this phenomenon in further detail, we 
introduced a Langevin equation (Q) with multiplicative 
and additive noise terms as a general model for a stochas- 
tic process of this type. We analyzed its stationary state 
theoretically and numerically, and found that this model 
actually reproduces the power law of moments. Further- 
more, by comparing the approximate and exact treat- 
ments of the effect of the additive noise, the usual ap- 
proximation of the additive noise by introducing a lower 
bound of the amplitude was justified. 

Although we restricted our study to the Langevin 
equation (yj) for the sake of precise argument, the power 
law itself is not sensitive to the details of the model em- 
ployed, and can appear robustly in many stochastic pro- 
cesses driven by multiplicative and weak additive noise. 
We demonstrated such robustness numerically for some 
different models, where the system is driven by noises 
which are not Gaussian-white. Furthermore, as some 
typical realizations of this type of power law, we dis- 
cussed the power law of moments in noisy on-off intermit- 
tency, and the power-law spatial correlation functions in 
the spatio-temporal chaotic regime of non-locally coupled 
systems, which gives some more insight into the power 
law of the spatial correlation function. 

As we already noted, the power law of moments seems 
to be a general phenomenon appearing in various systems 
over a wide range of parameters. This mechanism of 
generating a power law seems to be quite universal, and 
some of the power laws observed in the real world may 
belong to this class. 
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